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Downstream (molecular biology) 



Downstream refers to a relative position in DNA or RNA. Each 
strand of DNA or RNA has a 5" end and a 3' end, so named for 
the carbons on the deoxyribose ring. Relative to the position 
on the strand, downstream is the region towards the 3' end of 
the strand. Since DNA strands run in opposite directions, 
upstream on one strand is downstream on the other end. 

Another important meaning of the term "downstream", refers 
to the temporal and mechanistic order of cellular and 
molecular events. For example, in s igna l t ransduction , the 
second messenger acts downstream to activation of ceM 
membrane receptors. 
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Summary 

The transcriptional factor DREB/CBF (dehydration-responsive element/C-repeat-binding) specifically inter- 
acts with the dehydration-responsive element |DRE)/C-repeat (CRT) c/s-acting element (A/GCCGAC) and 
controls the expression of many stress-inducible genes in Arabidopsis. Transgenic plants overexpressing 
DREB1A showed activated expression of many stress-inducible genes and improved tolerance to not only 
drought salinity, and freezing but also growth retardation. We searched for downstream genes in trans- 
genic plants overexpressing DREB1A using the full-length cDNA microarray and Affymetrix GeneChip array. 
We confirmed candidate genes selected by array analyses using RNA gel blot and identified 38 genes as the 
DREB1A downstream genes, including 20 unreported new downstream genes. Many of the products of 
these genes were proteins known to function against stress and were probably responsible for the stress 
tolerance of the transgenic plants. The downstream genes also included genes for protein factors involved 
in further regulation of signal transduction and gene expression in response to stress. The identified 
genes were classified into direct downstream genes of DREB1A and the others based on their expression 
patterns in response to cold stress. We also searched for conserved sequences in the promoter regions 
of the direct downstream genes and found A/GCCGACNT in their promoter regions from -51 to -450 as 
a consensus DRE. The recombinant DREB1A protein bound to A/GCCGACNT more efficiently than to 
A/GCCGACNA/G/C. 

Keywords: abiotic stress, transcription factor, downstream genes, stress tolerance, A/GCCGACNT. 



Introduction 

Plants respond and adapt to environmental stresses such 
as drought, salinity, and cold with physiological and devel- 
opmental changes. These stresses induce many plant 
genes, which function in stress tolerance (for reviews, 
see Bray, 1997; Ingram and Battels, 1996; Shinozaki and 
Yamaguchi-Shinozaki, 1997, 2000). The products of these 
genes are thought to function not only in protecting cells 
from stress but also in the regulation of genes for signal 
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transduction in stress response (Shinozaki and Yamaguchi- 
Shinozaki, 2000). Thus, these gene products are classified 
into two groups. The first group includes proteins that 
probably function in stress tolerance: water channel pro- 
teins, the enzymes required for the biosynthesis of various 
osmoprotectants, late-embryogenesis-abundant (LEA) pro- 
teins that may protect macromolecules and membranes, 
enzymes for modifying membrane lipids, detoxification 
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enzymes, and so on, The second group contains protein 
factors involved in further regulation of signal transduction 
andgene^exf^ 



2000; Xasuga eiaL s 19199; Liu ef ai, 1998). Hbwever/over- 
expressipn'^ DREB1A protein also caused severe 

growth retardation under normal growth conditions. Use 

Mcon- 



Analyses of drought-inducible genes have indicated the minimizes negative effects on plant growth (Kasuga ef a/., 
^existence oiABA-independent. and.ABA^ 

transduction cascades' between the initial signal of water . stress-inducible genes of DBEB1A using RNA gel blot ana- 
deficit and the expression of specific genes (reviewed by tysis. We also identified further downstream stress-induci- 
Bray, 1997; Ingram and Bartels, 1996; Shinozaki and Yama- ble genes of DREB1A using 1300 full-length Arabidopsis 
guchi-Shihozaki, 20001 Several stress-inducible genes, cDNAs microarray (Seki et ai, 2001). Twelve stress-induci- 
such as rd29A and cor15a, are induced through this ABA- ble genes were identified as downstream genes of DREB1 A, 
independent pathway. A c/s-acting element has been iden- and six of them were novel. These DREB1A downstream 
tifled-in the' promoter region of rd29A and is responsible for genes contained the DRE or DRE-related core motifs in their 
drought-, high-salt-, and cold-inducib^ . . • • V priDmoterrejgions^ Recently; Fowlerand Thbmashow(2002) 

This .'sfiniiflnrA' (TACCdikCATS n.moH.tK. H*hwW«V; ia'+^a .a* a * /^r-'irww-n ....... 



in?/ 



exacting" elements, named C-repeat (CRT) and low-tem- 
perature-responsive 'element (LTRE), : bothfeoritaining an 
A/GCCGAC core motif,, regulate, cold-inducible genes 
(Baker ot at. 1994> Hann 1QQfi« Thnmnch^., iqqq\ 
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and DREB2A, have 'been isolated' by the yeast one-hybrid 
screening method (Liu et ai, 1998; Stockinger et ai, 1997). 
The deduced amino acid sequences of these proteins 
showed. significant sequence similarity in the conserved 



array. We confirmed candidate genes selected by array 
analyses using RNA gel blot analysis. We identified 38 
genes as the downstream genes of DREB1A, 20 of which 
were unreported new downstream genes. The products of 




CRT sequence in Arabidopsis. We isolated two cDNA clones downstream genes of DREB1A from the identified down- 




induced by coldstress/and expression of the DREB2 genes ment; and (in) existence of the DRE core motif in the 

.is induced by-drought and high salt stresses, DREB1/CBFs promoter region. 

.are thought to function in cold-responsive gene expression, ' " : ; ""' 
whereas DREB2s are involvedin drought-responsive gene 

expression: , Results 

Strip^g4oIerance to freezjng;stress was observed in the 
transgenic Arabidopsis plants overexpressing CBFT 

virus (CaMV) 35S promoter faa^i&Q^ 
reported that overexpression of the DREB1A [CBF3\ cDNA 

strong expression of downstream stress-inducible; genes 
:8™|^ | g he r ;j tble ranc e jta 

drpug.ht>/.Hlg^ a^^fre^inig stresses ffiilmdu 
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transcription factor 



under control of the CaMV 35S promoter in transgenic 
Arabidopsis plants <35S:DREBrA transgenic plants) gave 
rise to strong constitutive expression of the stress-inducible 
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drought and high salt stresses. Overexpression of DREB1 A 
also resulted in growth retardation under normal growth 
conditions. Previously, we identified 13 downstream cold- 
regulated genes in transgenic plants using RNA gel blot and 
microarray analyses (Kasuga et al, 1999; Seki et al., 2001; 
Taji et ai, 2002). To determine more downstream cold- 
reguiated .genes of DREB1A in Arabidopsis, we performed 
array analysis using the approximately 7000 RAFL cDNA 
microarray and the approximately 8000 Affymetrix Gene- 
Chip array. Preparation of the RAFL cDNA microarray con- 
taining the approximately 7000 Arabidopsis full-length 
cDNAs was described previously by Seki et al. (2001). A 
total of 6454 independent Arabidopsis genes were included 
in the RAFL cDNA microarray. The Affymetrix GeneChip 
array contains 7483 independent Arabidopsis genes. We 
compared genes of the RAFL cDNA microarray and the 
Affymetrix GeneChip, and found 1919 genes in common 
(Figure 1). In this study, we examined the expression pro- 
files of 12 018 independent genes, which accounted for 
almost half of the Arabidopsis genes. 

The RAFL cDNA microarray analysis was carried out as 
described before by Seki ef al. (2001, 2002). We used two 
independent transgenic Arabidopsis plants overexpressing 
DREB1A, 35S:DREB1Ad, and 35S:DREB1Af. Both transgenic 
plants showed growth retardation and stress tolerance to 
drought, high salt, and cold like 35S:DREB1Ab (Kasuga 
etaL, 1999). In both transgenic plants, the DREB1A gene 
was strongly expressed under control condition. mRNAs 
from the 35S:DREB1Ad and 35S:DREB1Af transgenic plants 
and wild-type (pBI121) plants were used for the preparation 
of Cy3- and Cy5-labe!ed cDNA probes, respectively. These 




Figure 1. Classification of genes on the RAFL cDNA microarray and the 
Affymetrix GeneChip array. 

A total of 6454 independent Arabidopsis genes were included in the 
approximately 7000 RAFL cONA microarray. The Affymetrix GeneChip 
had approximately 8000 DMA oligonucleotide probe sets, and the number 
of independent Arabidopsis genes was 7483. We compared genes of the 
RAFL cONA microarray and the Affymetrix GeneChip, and found 1919 genes 
in common. Therefore, we searched a total of 12 018 genes, which con- 
tained almost half of the Arabidopsis genes. Seventeen and 18 genes were 
identified as the downstream genes of DREB1A by the RAFL cDNA micro- 
array and the Affymetrix GeneChip, respectively, and 9 Identified genes 
were in common. We searched 13 434 independent Arabidopsis promoters 
that contain A/GCCGACNT and identified additional 12 genes as down- 
stream genes of DREB1A. The numbers of identified downstream genes are 
shown in the parentheses. 



cDNA probes were mixed and hybridized with the cDNA 
microarray. To assess the reproducibility of the microarray 
analysis, we repeated the experiment three times. Hybridi- 
zation of different microarrays with the same mRNA sam- 
ples indicated a good correlation. We focused on the genes 
with expression ratios greater than five times in at least 
either 35S:DREB1Ad or 35S:DREB1Af transgenic plants as 
compared with the control plants. Analysis with the RAFL 
cDNA microarray revealed an increase in the level of 19 
gene transcripts under the non-stressed normal growth 
condition in the transgenic plants. Analysis with the Affy- 
metrix GeneChip array was carried out as described by 
Goda ef al. (2002). The DREB1A downstream genes were 
analyzed using the 35S:DREB1 Ad and control plants. Signal 
log ratios were calculated for the duplicate transgenic plant 
samples compared with each of the two control samples, 
generating four comparisons. We adopted genes with sig- 
nal log ratios greater than two times in all four compar- 
isons. Analysis with the Affymetrix GeneChip array 
revealed an increase in the transcript level of 32 genes 
under the control condition in transgenic plants. 

We selected 19 and 32 genes as candidates for down- 
stream genes of the DREB1A protein using the RAFL cDNA 
microarray and the Affymetrix GeneChip array, respec- 
tively. Among these 51 candidate genes, 9 were upregu- 
lated genes found in common by both array analyses. 
Among the 42 genes, 8 were already reported as down- 
stream genes of DREB1A (Kasuga et ai, 1999; Seki et al, 
2001; Taji ef a/., 2002) and 34 candidate genes were unre- 
ported at that time and were further analyzed by RNA gel 
blot hybridization. Twenty-six genes showed increased 
mRNA accumulation under the control condition in both 
35S:DREB1 Ad and 35S:DREB1 Af transgenic plants as com- 
pared with that in wild-type plants (Figure 2a). We con- 
firmed 17 and 18 genes as the downstream genes using 
RNA gel blot analysis in the 19 and 32 candidate genes 
selected by the RAFL cDNA microarray and the Affymetrix 
GeneChip array, respectively. As 9 genes were found in 
common by both array analyses, 26 genes in total were 
identified as the downstream genes of DREB1A (Figure 1; 
Table 1). 

The circadian clock has been reported to control expres- 
sion of the DREB1A/CBF3 gene and some of the identified 
downstream genes of DREB1A (Harmer ef al, 2000). There- 
fore, we analyzed the effects of the circadian clock on the 
expression of DREB1A and its downstream genes in com- 
parison with the effect of cold stress. We analyzed the 
expression of DREB1 A and all identified downstream genes 
every 3 h for 1 day under a non-stressed control condition 
and found that all genes except AtGolS3 and At2g23120 
were induced in the daytime. The expression patterns of 
DREB1A and the two downstream genes rd29A and cor15A 
are shown in Figure 2(b). The level of expression of 
these genes in the non-stressed condition was less than 
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Figure 2. RNA gel-blot analysis of the DREB1A and its downstream genes. 

(a) Expression of DREB1A and its downstream genes in the 35S:DREB1A 
transgenic and wild-type plants. Each lane was loaded with 10 jig of total 
RNA prepared from the transgenic and wild-type (pBI121) plants. Lane 1 
contains RNA prepared from the wild-type plants, lane 2 contains RNA 
prepared from the 35S:0REB1 Ad transgenic plants, and lane 3 contains RNA 
prepared from the 35S:0RE81Af transgenic plants. 

(b) Analysis of effects of the circadian clock on the expression of DREB1A 
and its downstream genes. Each lane was loaded with 10 \ig of total RNA 
prepared from the wild-type plants. The expression of DREB1A and its 
downstream genes was analyzed every 3 h for 1 day underthe non-stressed 
control condition (lanes 1-6). Induction of DREB1A and its downstream 
genes by cold stress was analyzed after starting treatment at two different 
time points (lanes 9-20). 



one-tenth of the maximum induction level under cold 
stress. Expression of AtGolS3 was not detected under 
the non-stressed condition during the 24-h study period 
(Figure 2b). At2g23120was expressed even underthe non- 
stressed condition and the level of expression did not 
change for 24 h (Figure 2b). Then, we analyzed induction 
of all downstream genes by cold stress starting treatment at 
two different time points and found that induction of all 
DREB1A and its downstream genes showed the same 
expression pattern regardless of starting time (Figure 2b). 
These results indicate that DREB1A and its downstream 
genes are clock controlled but that these effects are not 
important for induction of these genes by cold stress. 

Cold regulation of the identified downstream genes of 
DREB1A 

We previously reported that expression of the DREB1A/ 
CBF3 gene was regulated by cold stress. The expression 
level of the DREB1A/CBF3 transcript increased within 1 h 
after exposure to cold stress (4°C) and reached its max- 

© Blackwell Publishing Ltd, The Plant Journal, (2004), 38, 982-993 



imum at 2 h (Liu ef al., 1998; Figures 2b and 3). We verified 
whether the identified new downstream genes of DREB1A 
were also regulated by cold stress using RNA gel blot 
hybridization. As shown in Figure 3 f all 26 identified genes 
were induced by cold stress within 24 h. We compared 
the expression patterns of these downstream genes in 
response to cold stress and classified them into three 
groups: Class l-Class III. 

The Class I genes showed markedly increased expression 
within 5 h and reached a maximum at 10 or 24 h after 
exposure to cold stress. These Class I genes are indicated 
in blue in Figure 3. Many reported DREB1A downstream 
genes, such as rd29A, corWa, corWb, kin1, kin2, erdW, 
rd17, and AtGolS3 belong to this class. The expression of 
the Class II genes was induced within 2 h and reached a 
maximum at 5 h after exposure to cold stress. Only one 
gene belongs to this class and is shown in orange in 
Figure 3. The Class III genes were induced gradually up 
to 24 h. Eight genes belong to this class and are shown in 
green in Figure 3. The expression of the DREB1A gene 
during cold stress preceded that of the Class I genes, which 
have the DRE core motif in their promoter (except for 
At4g12470). These results indicated that the Class I genes 
may be direct downstream genes of DREB1A. On the other 
hand, the timing of the expression of the Class II and Class 
III genes is different from that of the Class I genes. These 
genes probably were regulated by DREB1A through other 
factors, such as DREBlA-induced transcription factors in 
Arabidopsis, or regulated by both DREB1 A and other tran- 
scription factors. 

Conserved sequence in the promoter region of the 
downstream genes of DREB1A confirmed by 
RNA gel blot analysis 

In our previous study, the DREB1AICBF3 gene encoding 
DRE-binding protein was isolated by the yeast one-hybrid 
screening method using Arabidopsis cDNA libraries (Liu 
et ai, 1998). The DREB1A protein specifically bound to six 
nucleotides (A/GCCGAC) of the DRE core motif (Sakuma 
ef a/., 2002). We searched for conserved sequences in the 
promoter regions of the 26 downstream genes of DREB1A 
confirmed by RNA gel blot analysis using the meme 
(Bailey and Elkan, 1995) and SEQUENCE LOGO programs • 
(Schneider and Stephens, 1990), The 1000-bp promoter 
regions upstream of putative transcriptional initiation sites 
that are the 5' ends of the full-length cDNA clones were used 
for this search (Seki ef a/., 2002). We found two highly 
conserved sequences, A/GCCGACNT and ACGTGG/T, that 
included the DRE core and ABRE (ABA response element) 
core motifs in these 1000-bp promoter regions (Figure 4a). 
Among the 26 downstream genes, 18 (69%) and 17 (65%) 
contained the A/GCCGACNT and ACGTGG/T sequences in 
their 1000-bp promoter regions, respectively. In addition, 
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Table 1 List of genes that are upregulated by 0REB1A overexpression and cold stress 

RAFL clone AGI Gene name Description Methods Publication 



Transcriptional factor 



RAFL04-15-K19 


At 1g 27730 


S 77 ZAT 10 


C 2 H 2 zinc finger DNA-binding protein 


Promoter 


New 


RAFL21-19-J24 


At5g04340 




C 2 H 2 zinc finger DNA-binding protein 


Promoter 


New 




At1g46768 


RAP? 1 


Mrt uumain uixjM-utfiuiny pruiein 


r romoicr 


A 

*+ 


PI metabolism 












RAFL09-43-P14 


At5g58670 


AtPId 


Phosphatidylinositol-specific phospholipase 


8000 


New 


RNA-binding protein 










RAFL04-16-B18 


At2g21660 


AtGRP7 


Glycine-rich RNA-binding protein 


8000 


New 


Transport protein 












RAFL09-48-J19 


At4g35300 






Prnmntpr 
i lUlllUlGf 


New 


Hoc a t i i ra c a 












RAFL05-09-K18 


At2g31360 


ADS2 


Delta 9 desaturase 


8000 


New 


ivieiaooHsrn 












RAFL04-16-K22 


At 1g 093 50 


All UUIOJ 




7onn flooo 

/uuu, ouuu 


3 4 


RAFL19-51-F16 


At4g23600 




Tyrosine transaminase 


8000 


4 


RAFL09-23-E17 


At4g33070 


Pdcl 


Pyruvate decarboxylase-! 


Promoter 


New 


Protease inhibitor 












RAFL05-08-E12 


At2g02100 




Gamma-thionins family 


7000 


New 


RAFL05-18-K08 


At2g43510 




Tryosin inhibitor 


8000 


New 


RAFL1MVM23 


At4g 12470 




Protease inhibitor/seed storage/LTP family 


7000, 8000 


4 


RAFL06-76-D14 


At4g 12480 




Protease inhibitor/seed storage/LTP family 


8000 


New 


LEA protein 












o api nK.n an 


At 1 nrt1 A~ll\ 

At igui*»/u 




Late embryogenesis abundant protein 


7000, 8000, promoter 


4 


bapi nn ?i Pit 

rlArLUD-Z It Ij 


At 1 <i 1 cocn 
At ig I bBDU 




Laic tjniuryugeiitibib duunoani proiein 


7nnn 
/uuu 


New 


DA CI f\A On MHO 


At 1 m tf\AAI\ 

At igzU44U 


tQ I / cur **/ 


i/Bnyann 


i^nn 7nnn ftnnn 

IOUU, /uuu, ouuu 


1 9 , 


nArLUO-Uo-rl / 


Ati n nc\A en 

At igzo4bu 


ciu fu 


L/enyanii 


1 **nn 7nnn Rnnn 
uuu, /uuu, ouuu 


I, ' 


R api np on Mir 


AtOn/icin 
Atiiga/o JU 


cor 7 5b 


Laic ciiiui yuytJiicaia dUUnUdlil piuicin 


7nnn ftnnn nrnmntcr 

/uuu, ouuu, pruiiiuici 


4 


DACI ntM AHP 


AtOz-i/ioc/i n 
AI^g4za4U 


cor 1 5a 


1 fitp pmhrunnpnpik ahimrlant nrntpin 
Laic ciiiui yuyciicoio duuiiuaiu piuiciti 


1 JUU, /uuu 


1 2 


nMrLUs* 1 3*nU3 


Mioy juj /u 




ucnyui til 


ftnnn nrnmnfpr 

ouuu, prumuicr 


A 
H 


RAFL14-61-K05 


At4g15910 




1 ato omkn/nnonocic ahnnWant nrntoin 
Ldlo clllUi yuyclicala dUUIiUdlu \j\ Olcln 


P rn m aJp r 

r rornoicr 


A 


KIN protein 












RAFL06-08-N16 


At5g15960 


kinl 


Late embryogenesis abundant protein 


1300, 7000 


1,2 


RAFL04-17-B12 


At5g 15970 


kin2 cor6.6 


Late embryogenesis abundant protein 


1300, 7000 


1,2,- 


HydrophiNc protein 












RAFL04-17-F01 


At5g52310 


rd29A Iti78 


HvrtrnnhNir nrntpin 

1 1 y u i w yj i ii t lit piuiwiii 


Hon 700D 8000 Drnmoter 


1 2 - 


Senescence-related 


protein 










RAFL08-19-H17 


At2g 17840 


erd7 


Unknown protein 


Promoter 


4 


Unknown protein 












RAFL05-1 8-020 


At1g29390 




Unknown protein 


7000 


New 


RAFL04O2-F24 


At1g51090 




Proline-rich protein 


Promoter 


New 


RAFL06-15-K18 


At1g54410 




Unknown protein 


7000, 8000 


New 


RAFL04-10-D13 


At2g23120 




Unknown protein 


7000 


4 


RAFL05-14-K23 


At2g28900 




Unknown protein 


8000 


4 


RAFL02-02-B06 


At2g33830 




Unknown protein 


7000, 8000 


New 


RAFL09-65-H03 


At4g 14000 




Unknown protein 


Promoter 


New 


RAFL07-07-J02 


At4cj24960 


AtHVA22d 


HVA22 family 


Promoter 


New 


RAFL06-07-E01 


At4g30650 




Unknown protein 


7000 


New 


RAFL03-01-H06 


At4g38580 


ATFP6 


Farnesylated protein 


Promoter 


New 


RAFL08-11-P07 


At5g174€0 




Unknown protein 


Promoter 


New 


RAFL0<M8-807 


At5g62350 




Unknown protein 


1300, 8000 


2, 4 



Publication indicates whether the gene has not been reported as a DREB1 A downstream gene before (New) or has already been reported as 
a DREB1A downstream genes. 1, Kasuga et al (1999); 2, Seki et al. (2001); 3, Taji et al. (2002); and 4, Fowler and Thomashow (2002). 
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Figure 3. Expression of the DREB1 A downstream genes in response to cold 
stress. 

Each lane was loaded with 10 \ig of total RNA prepared from Arabidopsis 
plants that had been transferred to and grown for the indicated time at 4 B C. 
The intensity of each band quantified by densitometry shows the levels of 
expression of the downstream genes. Relative signal values at 2 and 10 h for 
each gene ere indicated by red squares. We compared the expression of 
these downstream genes in response to cold stress and classified them into 
Class (-Class III. The Class I genes showed markedly increased expression 
within 5 h with a maximum at 10 or 24 h, after exposure to cold stress. These 
Class I genes are indicated by gene names in blue. Expression of the Class II 
genes was induced within 2 h and reached a maximum at 5 h after exposure 
to cold stress. One gene belongs to this class and is shown by gene name in 
orange. The Class ill genes were induced gradually during 24 h. Eight genes 
belong to this class and are shown by gene names in green. 
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among the 17 Class I genes, 16 (94%) and 14 (82%) con- 
tained the A/GCCGACNT and ACGTGG/T sequences in their 
promoters, respectively. 

We estimated the frequency of these conserved 
sequences in the 1000-bp promoter region in 13 434 inde- 
pendent genes of Arabidopsis and compared it with that in 
the 1000-bp promoter region of the 26 downstream genes. 
First, we examined the base contents of the 1000-bp 
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Figure 4, Analysis of DRE and ABRE in the promoter regions of theDREBIA 
downstream genes or the 13 434 genes in Arabidopsis. 

(a) Sequence logo for the DRE and ABRE core sequences in the promoter 
regions of the 26 DREB1A downstream genes. DRE and ABRE core 
sequences in the promoter regions of the 26 DREB1A downstream genes 
were used to create the sequence logo. The relative frequencies of bases are 
given by their relative heights. Statistical significance of the 13th nucleotide 
T is determined by x 2 test [P < 0.001). 

(b) The base composition of the 13 434 Arabidopsis promoter sequences. 

(c) The base composition of 13 434 Arabidopsis promoter sequences in 
every 50 bp of the 1000-bp promoter sequences. The GC content of Arabi- 
dopsis promoter is 33.4%. 

(d) The number of DRE and ABRE in every 50 bp of the 1000-bp promoter 
sequences of the 13 434 genes. The expected values of DRE and ABRE are 
indicated by line charts. Statistical significance was determined by x 2 test. 
There are significant differences in frequency of the DRE core and ABRE core 
motifs in -450 to -51 promoter regions between the found value and the 
expected value (P < 0.001), respectively. 

(e) The numbers of the genes having the DRE core and ABRE core motifs in 
their promoter regions from -450 to -51, are indicated by colored circular 
charts. The numbers of genes with the A/GCCGAC, A/GCCGACNT, and 
ACGTGG/T sequences in their promoters are indicated by sky blue, blue, 
and orange, respectively. DRE and ABRE core elements were found more 
frequently than their expected number in the promoter regions from -51 to 
-450 of the 13 434 Arabidopsis genes. 

Predicted promoter regions of 26 DREB1A downstream genes were com- 
pared by the meme and sequence logo programs. The 1000-bp pro- 
moter regions upstream of putative transcriptional initiation sites that are 
the 5' ends of the full-length cDNA clones were used for this analysis. 



promoter sequences in the 13 434 genes and estimated 
the expected value of these elements. Then, we compared 
the expected value with the found value of these two 
elements in the 13 434 promoters. The promoter sequences 
in the database were AT rich and A:T:G:C was 
4 533 367 : 4 407 824 : 2 209 509 : 2 283 300, respectively. 
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The GC content of the Arabidopsis 13 434 promoters was 
33.4% (Figure 4b). The expected values of the A/GCCGA- 
CNT and ACGTGG/T sequences were 1163 and 6954, 
respectively. On the other hand, the found values of 
A/GCCGACNT and ACGTGG/T were 1703 and 9610, respec- 
tively. There was no significant difference between the 
expected value and the found value of the frequency of 
A/GCCGACNTinthe13 434 independent promoters^ < 0.05; 
X 2 test). However, the frequency of ACGTGG/T was signifi- 
cantly different between these two values in the 13 434 
promoters (P < 0.001; x 2 test). The found value was higher 
than the expected value of the frequency of A/GCCGACNT. 

Then, we divided the 1000-bp promoter sequences into 
every 50 bp, and searched for the base contents of each 
50 bp (Figure 4c). The expected numbers of A/GCCGAC and 
ACGTGG/T were calculated in the base contents of every 
50 bp of the 1000-bp promoter sequences (Figure 4d). We 
found a hotspot in both DRE core and ABRE core motifs in 
the promoter region from -51 to -450 bp upstream of 
putative transcriptional initiation sites of the 13 434 Arabi- 
dopsis genes (Figure 4d). The found value of the two motifs 
was significantly higher than the expected value in the 
promoter region from -51 to -450 bp upstream of the 
putative transcriptional initiation sites (P < 0.001; x 2 test). 

We searched for A/GCCGAC, A/GCCGACNT, and 
ACGTGG/T in the -51 to -450 promoter region of the 26 
downstream genes of DREB1A and the 13 434 independent 
genes. We found A/GCCGAC, A/GCCGACNT, and ACGTGG/ 
T located in the -51 to -450 promoter regions of 18 (69%), 
16 (62%), and 14 (54%) of the downstream genes, respec- 
tively (Figure 4e). On the other hand, we found A/GCCGAC, 
A/GCCGACNT, and ACGTGG/T located in the -51 to -450 
promoter region of 2270 (16.9%), 531 (4.0%), and 3674 
(27.3%) genes of the 13 434 Arabidopsis genes, respectively 
(Figure 4e). The frequency of these sequences was signi- 
ficantly different between the 26 identified genes and the 
13 434 independent genes (P < 0.001; x z test). These 
results indicate that most of the downstream genes have 
DRE and ABRE in the -51 to -450 promoter region. In 
addition, the frequency of the A/GCCGACNT sequence in 
the promoter region of the 26 downstream genes was very 
high (69%), especially in that of the 17 Class I genes (94%). 
The A/GCCGACNT sequence may be important for the 
direct downstream genes of DREB1A. 

Competitive DNA binding assay of DREB1A 

We performed a gel mobility shift assay with DNA frag- 
ments containing the ACCGACAT or ACCGACAA/C/G 
sequences as competitors to analyze binding affinity of 
DREB1A. The extent of binding of DREB1A was reduced 
more effectively by the addition of the unlabeled DNA 
fragment containing ACCGACAT than by the addition of 
the unlabeled DNA fragment containing ACCGACAA/C/G 
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Figure 5. Characterization of the DNA-binding affinities of the recombinant 
DREB1A proteins. 

Sequences of the 75-bp fragment of the wild-type rd29A promoter (ACCGA- 
CAT) and the mutated fragments (ACCGACAG, ACCGACAA, and ACCGA- 
CAC) used as competitors. The ACCGACAT fragment was also used as a 
probe. Gel mobility shift assay of DREB1A showing a different binding 
affinity between the wild-type and mutated competitors. GST (glutathion 
S-transf erase) or recombinant protein solutions were pre-incubated with or 
without competitors for 5 min at 25*C Then, the M P-labeled probe was 
added and the mixture was incubated for 30 min at 25°C As competitors, 
10-fold (xlO), 100-fold (xlOO), or 1000-fold (x1OO0) excess amounts of the 
unlabeled fragments were used. 



(Figure 5). These observations confirmed the results 
obtained from the analysis of the promoter region of the 
DREB1A downstream genes (Figure 4). DREB1A binds pre- 
ferentially to ACCGACAT rather than to ACCGACAA/C/G. 

Downstream genes of the DREB1A protein identified by 
using Arabidopsis promoter sequence data and RNA gei 
blot analysis 

As the DRE core motif was found frequently in the promoter 
region of downstream genes of DREB1A, we searched for 
Arabidopsis genes that contained multi-A/GCCGACNT and 
-ACGTGG/T to identify other downstream genes. We found 
15 genes having multiple elements of both A/GCCGACNT 
and ACGTGG/T in their promoters. We selected these genes 
as candidate DREB1A downstream genes. Five of the 15 
genes have already been identified as downstream genes. 
Therefore, we analyzed the expression of the 10 genes in 
the 35S:DREB1Ad and 35S:DREB1Af plants using RNA gel 
blot hybridization. Among them, four genes were over- 
expressed in both transgenic plants under a control condi- 
tion as compared with the wild-type plants (Figure 6; 
Table 2). 

To analyze the downstream genes of DREB1A compre- 
hensively, we searched further for A/GCCGACNT in the 
promoter region of cold-regulated genes that were identi- 
fied using array analyses (Chen et al., 2002; Fowler and 
Thomashow, 2002; Seki et al., 2002). We selected 30 genes 
as candidate DREB1A downstream genes and analyzed 
expression of these genes in the transgenic plants using 
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Figure 6. RNA gel-blot analysis of the downstream genes of DREB1A iden- 
tified by promoter analysis and their expression in response to cold stress. 

(a) Expression of the downstream genes of DREB1A in the 35S:DREB1A 
transgenic and wild-type plants. Each lane was loaded with 10 of total 
RNA prepared from the transgenic and wild-type (pBI121) plants. Lane 1 
contains RNA prepared from the wild-type plants, lane 2 contains RNA 
prepared from the 35S:DREB1Ad transgenic plants, and lane 3 contains RNA 
prepared from the 35S:DREB1Af transgenic plants. 

(b) Each lane was loaded with 10 ug of total RNA prepared from Arabidopsis 
plants that had been transferred to and grown for the indicated lime at 4°C. 
The intensity of each band quantified by densitometry shows the levels of 
expression of the downstream genes. Relative signal values at 2 and 10 h for 
each gene are indicated by red squares. We compared the expression of 
these downstream genes in response to cold stress and classified them into 
Class l-Class III. The Class I genes showed markedly increased expression 
within 5 h with a maximum at 10 or 24 h, after exposure to cold stress. These 
Class I genes are indicated by gene names in blue. Expression of the Class II 
genes was induced within 2 h and reached a maximum at 5 h after exposure 
to cold stress. One gene belongs to this Class and is shown by gene names 
in orange. The Class III genes were induced gradually during the 24 h. One 
gene belongs to this class and is shown by gene names in green. 



RNA gel blot hybridization. Eight genes were overex- 
pressed in transgenic plants under the control condition 
(Figure 6a; Table 2). We could not detect overexpression of 
the other 22 genes either in the transgenic plants or the 
wild-type plants exposed to cold stress (data not shown). 
Therefore, we identified these additional 12 genes as the 
downstream genes of DRE81A using the promoter 
sequence data of Arabidopsis genes and RNA gel blot 
analysis. Twelve additional downstream genes of DREB1A 
were classified into three groups based on their expression 
patterns in response to cold stress. Among the 12 genes, 10 
belong to Class I, 1 belongs to Class II, and the other 1 
belongs to Class III (Figure 6b). A total of 38 genes were 
identified as the downstream genes of DREB1A and these 
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genes included 20 unreported new downstream genes 
(Table 1). 



Characterization of the downstream genes of DREB1A 

We indicated descriptions of the identified genes based 
on the results of blast (Altschul et a/., 1990) and pfam 
.(Bateman et a/., 1999) searches in Table 1. These gene 
products are transcription factors, phospholipase C, RNA- 
binding protein, sugar transport protein, desaturase, car- 
bohydrate metabolism-related proteins, LEA proteins, KIN 
(cold-inducible) proteins, osmoprotectant biosynthesis- 
protein, protease inhibitors, and so on. Twenty-three 
identified genes have been reported as the DREB1A down- 
stream genes previously and 20 genes were newly identi- 
fied as the downstream genes in this study. Transcriptional 
factors contain two kinds of C 2 H 2 -type zinc-finger DNA- 
binding proteins. An RNA-binding protein has an RNA- 
recognition motif containing a glycine-rich region. We have 
identified genes for many kinds of LEA proteins and pro- 
tease inhibitors as the downstream genes. The functions of 
10 identified genes remain unknown (Table 1). 

Discussion 

We searched for downstream stress-inducible genes of 
DREB 1A in 35S:DREB1A plants using the RAFL cDNA micro- 
array and the Affymetrix GeneChip array. We selected 19 
candidate genes using the RAFL cDNA microarray and 
confirmed 17 of them as the downstream genes of DREB1A 
by RNA gel blot analysis. In addition, we selected 32 can- 
didate genes using the Affymetrix GeneChip array, and 
confirmed 18 of them as the downstream genes of DREB1A 
by RNA gel blot analysis. The RAFL cDNA microarray con- 
tains cDNAs isolated from cDNA libraries prepared from 
Arabidopsis plants exposed to stress such as drought, high 
salt, and cold stresses. Therefore, most of the stress-indu- 
cible genes may be included in the cDNA microarray. To 
select the downstream stress-inducible genes showing low 
expression or that are not included in the microarrays, we 
searched the Arabidopsis promoter sequences and identi- 
fied 12 additional genes as the downstream genes of 
DREB1A. A total of 38 genes were identified as DREB1A 
downstream genes, and these genes included 20 new down- - 
stream genes. All 38 downstream genes were also upregu- 
lated by cold stress in Arabidopsis plants (Figures 3 and 6; 
Kasuga et al, 1999; Seki et a/., 2001). Most of the DREB1A 
downstream genes may be included in these 38 genes. 

The downstream genes of DREB1A have already been 
detected using microarrays by Seki et al. (2001) and Fowler 
and Thomashow (2002), previously. Fowler and Thoma- 
show used the Affymetrix GeneChip and indicated that 41 
genes were upregulated in all the CBF1/DREB1B-, CBF2/ 
DREB1C-, and CBF3/DREB1A-overexpressing transgenic 
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Table 2 List of genes that are found by promoter analysis 



RAFL clone 


AGI 


Gene name 


DRE 


Position 


ABRE 


Positio 


DREB1A downstream genes, which have multiple-DRE and ABRE elements 
RAFL04-12-F24 At1g51090 GTGAC G CCG ACATC AA . 

AATGAGCCGACATGGA 


-202 
-341 


AAATGACGTGTATGAA 


-141 


RAFL04-17-FQV 


At5g52310 


rd29A Iti78 


C ATG G A C CGACTACTA 
ATCATACCGACATCAG 
ATACTACCGACATGAG 
ATCAAGCCGACACAGA 


-275 
-225 
-168 
-131 


TTCATACGTGTCCCTT 


-65 


RAFLG5-03-A05* 


At2g42540 


cor15a 


TGTTG GCCGACATACA 
TCATGGCCGACCTGCT 
AATAAACCGACAAGGT 


-360 
-183 
-417 


ATTACACGTGGCCTGA 
AGGCCACGTGTAATCA 
CCTTCACGTGTATTTA 
TTAATACGTGTAATTT 


-131 
-128 
-304 
-70 


RAFL05-17-B13* 


At1g01470 




ATTC C ACCG ACGTGCA 
CATCGACCGACTTCAT 
AATGGACCGACCATGT 


-414 
-66 
-384 


ATGTCACGTGTTGATG 
GTTATACGTGTCTTCT 
TGTGTACGTGTGAGGG 


-646 
-147 
-103 


RAFL06-13-P08' 


At3g24190 




GTGGCGCCGACGTAGC 
TCGTTGCCGACGTAAT 


— OUo 

-263 


TTnrrACGTGGACCCT 


-100 


RAFL07-07-J02 


At4g 24960 


AtHVA22d 


CTCTCACCGACCGACG 
CACCGACCGACGTCTT 
CTGATGCCGACATACA 
ATG GC ACCG ACCTTAA 


-335 

— OO 1 

-184 
-297 


CCCACACGTGGCGAGA 
TCGCCACGTGTGGGAT 


-259 
-262 


RAFL08-11-P07 


At5g 17460 




GCCACGCCGACATAGT 
AACAG GCC GACATAAT 


-325 
-139 


CTGAAACGTGTCTATC 
TAAACACGTGGCTTAA 
AAGCCACGTGTTTAGT 
ACCCTACGTGGAAACA 


-204 
-233 
-236 
-81 


RAFL09-19-H05* 


At3g50970 


Xero2 


ACTACACCGACGTCTT 
TGTTGGCCGACATCGT 
TAGCGACCGACGTAAT 


-237 
-162 
-254 


TCAATACGTGTTGGCC 


-154 


RAFL21-19-J24 


At5g04340 




AAGTAGCCGACTTAAT 
TCTTAGCCGACTTCCA 


-412 
-250 


CAAACACGTGTACCAA 


-234 


Non-DREB1A downstream genes, which have multiple DRE and ABRE elements 
RAFL05-12-H13 At4g01020 AAGACACCGACATATA 

AATTCACCGACTTCAT 
GATGAACCGACAAATC 


-75 
-404 
-112 


TTGGAACGTGTAAATC 
TAAGCACGTGTGTATA 
AGATTACGTGTCTTAG 


-395 
-196 
-136 


RAFL06-80-P19 






GTTTTACCGACTTGCT 
G ATTTGCCG ACTTG G G 
GATTTGCCGACTTGGG 


-51 
-174 
-87 


AAATCACGTGTCCGAA 


-81 


RAFLQ8-12-K03 


At2g 15320 




GCAAAGCCGACATCCA 
GCAAAGCCGACCTCCA 


-296 
-81 


TACTAACGTGTTATTG 
AAAGTACGTGGTGGAG 


-455 
-65 


RAFL14-26-J19 


At3g46640 




CAAACGCCGACGGTGA 
TTACTACCGACTTAGC 
TTCAAACCGACCTAAC 


-852 
-293 
-125 


TTCCAACGTGGACGAG 
ACTCCACGTGGCTCCA 
TGATGACGTGTTGAAA 


-76 
-240 
-154 


RAFL15-18-E07 


At2g39810 


HOS1 


GATTTGCCGACTTGGG 
GATTTGCCGACTTGGG 
GTTTTACCGACTTGCT 
AAATAGCCGACCTGAA 


-234 
-147 
-270 
-336 


AAATPACGTGTCCGAA 


-240 


RAFL16-35-J05 


At3g 50960 




TAGCGACCGACGTAAT 
TGTTGGCCGACATCGT 
ACTACACCGACGTCTT 


-166 
-258 
-183 


TCAATACGTGTTGGCC 


-266 


DREB1A downstream genes, which are selected by low-temperature and DRE elements 
RAFL03-01-H06 At4g38580 ATFP6 GTCACACCGACGTACC -480 

G AC GTACC G AC CTTG A -488 






RAFL04-15-K19 


At1g27730 




TTATAGCCGACCTCTT 


-285 


CCGAAACGTGTACCAT 
CACACACGTGTACTAG 


-400 
-240 
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Table 2 continued 

RAFL clone AG! Gene name DRE Position ABRE Position 



RAFL08-19-H17 


At2g 17840 


erd7 


AAGCGACCGACCGACA 
GACCGACCGACA 1 GAG 


-210 
-ZQb 


TTTTAACGTGGCAATC 
n. n a a fx a r f5Tf5 fi at c TV 

OuMMOMUO 1 UUn lull 


-989 


RAFL09-23-E17 


At4g33070 


Pdd 


TTTAGACCGACATAAA 


-187 


GTTATACGTGGGAACA 
GATCTACGTGTATCTT 


-339 
-694 


RAFL09-48-J19 


At4g35300 




ATTATGCCGACATTAA 


-389 






RAFL09-65-H03 


At4g 14000 




CACAGACCGACTTTAA 
CTCGTACCGACCGGTT 
TGGAAGCCGACTAAAA 
CGGAAGCCGACCAAAG 
TTTCTACCGACTGAAT 


-994 
-62 
-658 
-571 
-398 






RAFL14-61-K05 


At4g15910 




ATCTCACCGACCTCTT 


-388 


CTAGAACGTGTCACTT 


-413 


RAP2.1 


At1g46768 




CTATGGCCGACATAAA 
AGGAGACCGACTGAAA 


-44 
-149 







•The genes with an asterisk mark are identified in array analysis. 



plants and 30 of them were cold-upregulated genes. Only 
11 of the 41 genes identified by Fowler and Thomashow 
were identical to those of the 18 downstream genes that we 
identified by Affymetrix GeneChip array analysis. This dis- 
crepancy may have been caused by the difference in the 
annotation table used for Affymetrix GeneChip array ana- 
lysis. The annotation tables of the Affymetrix GeneChip 
array have been frequently updated. We used the TAIR (The 
Arabidopsis Information Resource) annotation table 
updated in December 2002. In addition, we used the Arabi- 
dopsis thaliana ecotype Columbia, but they used A thali- 
ana ecotype Wassilewskija-2. Different stress-inducible 
downstream genes of DREB1A may function in these dif- 
ferent ecotypes. However, it may be difficult to determine 
the downstream genes accurately using only array analysis 
with repetition experiments. Therefore, we confirmed the 
downstream genes of DREB1A using RNA gel blot analysis. 
Only 18 genes were confirmed as DREB1A downstream 
genes by RNA gel blot analysis among 32 genes identified 
by the Affymetrix GeneChip. 

The direct downstream genes for DREB1A/CBF3 satisfied 
the following three criteria after RNA gel blot analysis: (i) 
increased level of transcripts in the 35S:DREB1A transgenic 
plants under a control condition; (ii) significantly increased 
transcripts between 2 and 10 h after cold treatment; and (iii) 
existence of the DRE core motif in the promoter region. The 
expression of the direct downstream gene was increased 
markedly after 2 h of cold treatment as the expression of the 
DREB1A/CBF3 transcript reaches its maximum at 2 h after a 
temperature shift from 22 to 4°C (Liu et aL, 1998; Figure 2b). 
In the present study, expression levels of the Class I genes 
were increased markedly within 5 h and reached their 
maximum at 10 or 24 h after exposure to cold stress 
(Figure 3). All the Class I genes except At4g12470 contained 
a conserved sequence (A/GCCGACNT) as a DRE core 



sequence in their -51 to -450 promoter regions, while 
most of the Class II and Class III downstream genes did 
not contain this sequence in their promoter regions. These 
results indicate that the Class I genes excluding At4g12470 
are the direct downstream genes of DREB1A. 

To investigate the importance of A/GCCGACNT as a DRE 
consensus sequence, we searched for the genes that have 
multiple A/GCCGACNT and ACGTGG/T in the -51 to -450 
promoter region. We selected 15 candidate genes in Ara- 
bidopsis (Table 2). Nine of them were upregulated under 
the control condition in the 35S:DREB1A transgenic plants, 
and were shown to be upregulated by cold stress using 
RNA gel blot analysis (Figure 6). However, the other six 
genes were not overexpressed in the transgenic plants 
(data not shown). In addition, we performed another search 
to examine the role of A/GCCGACNT on cold-inducible 
gene expression regulated by DREB1A. We searched for 
the genes containing A/GCCGACNT in the -51 to -450 
promoter region that were upregulated by cold stress. 
The cold-inducible genes were selected using array analy- 
sis data demonstrated by Seki et aL (2001, 2002), Fowler 
and Thomashow (2002), and Chen et aL (2002). We selected 
30 candidate genes in Arabidopsis, and analyzed the 
expression of these genes under a cold stress condition 
using RNA gel blot analysis. Eight of them were overex- . 
pressed in the 35S:DREB1A plants and upregulated by cold 
stress using RNA gel blot analysis (Table 2; Figure 6). How- 
ever, we could not detect the overexpression of 22 genes 
either in the 35S:DREB1A plants or in the wild-type plants 
exposed to cold stress (data not shown). Furthermore, we 
found that 531 genes contain A/GCCGACNT in their pro- 
moter regions from -450 to -51 in the 13 434 Arabidopsis 
genes. Many of these genes were not influenced by 
DREB1A, although they had A/GCCGACNT. These results 
indicate that A/GCCGACNT in -450 to -51 promoter region 
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is an efficient c/s-elements for the expression of genes 
regulated by DREB1A but another element, such as the 
other specific sequences or some secondary structure of 
the promoters, in the promoter region may be necessary for 
the expression of these genes. 

We investigated the functions of the identified genes 
(Table I), and estimated how overexpression of DREB1A 
in the transgenic plants improved stress tolerance. The 
identified downstream genes of DREB1A were classified 
into two groups. The first group includes proteins that 
probably function in stress tolerance such as LEA proteins, 
antifreeze proteins, hydrophilic protein, RNA-binding pro- 
tein, an enzyme required for biosynthesis of sugars (galac- 
tinol synthase; Taji et al., 2002), and protease inhibitors. 
These gene products probably function to increase the 
tolerance to drought, high salt, and freezing stresses in 
transgenic Arabidopsis. 

The second group contains protein factors involved in 
further regulation of signal transduction and gene expres- 
sion that probably functions in response to stress. They are 
transcription factors (C 2 H 2 zinc finger DNA-binding proteins 
and ERF/AP2-type DNA-binding protein) and enzymes in 
phospholipid metabolism (phospholipase C). The genes 
that encode the transcriptional factors containing C 2 H 2 -type 
zinc-finger motif constitute a large family in higher plants 
(Takatsuji, 1999). The transcription factors STZ and 
At5g04340 are such DREB1A downstream genes. STZ has 
been shown to repress the trans-activation of genes 
through an essential DLN-box/EAR-like repression motif 
in its C-terminal region (Ohta et al, 2001). The At5g04340 
protein also has this repression motif in its C-terminal 
region. Recently, we analyzed transgenic Arabidopsis 
plants overexpressing STZ using a cDNA microarray and 
found that many genes related to photosynthesis and car- 
bohydrate metabolism were suppressed in these plants 
(unpublished results by Sakamoto eta/.). These results 
suggest that STZ is involved in the mechanism of growth 
retardation of the 35S:DREB1A transgenic plants. Stress- 
related signaling has been reported to result in the repres- 
sion of genes associated with plant growth and develop- 
ment (Logemann eta/., 1995). Thus, DREB1A controls not 
only activation of genes involved in stress tolerance but 
also repression of genes involved in plant growth and 
development, such as photosynthesis-related genes. The 
wide range of functions of the downstream genes of 
DREB1A indicates that DREB1A regulates a complex gene 
expression network in response to cold stress. 

Experimental procedures 

Plant materials and cold stress treatments 

Arabidopsis thaliona (ecotype Columbia) and the transgenic 
35S:DREB1A plants {Kasuga et al., 1999) were grown under 



conditions of 16 h light/8 h darkness at 22°C on germination 
medium agar plates for 3 weeks. Cold stress treatment was per- 
formed as described before by Yamaguchi-Shinozaki and Shino- 
zaki (1994). 



Microarray analysis 

Total RNA was prepared using TRIZOL Reagent (Life Technologies, 
Rockville, MD, USA) from whole Arabidopsis plants, and mRNA 
was prepared using the mRNA isolation kit (Miltenyi Biotec, CA, 
USA), as described by Seki et al. (2001). For the RAFL cDNA 
microarray analysis, each mRNA sample was reverse transcribed 
in the presence of Cy3-dUTP or Cy5-dUTP (Amersham Pharmacia, 
Piscataway, NJ, USA). Hybridization and scanning analysis were 
carried out as described before (Seki et al., 2002). For the analysis 
using the Affymetrix GeneChip array, biotin-labeled cRNA probe 
sets were generated from the double-stranded cDNAs using the 
BioArray High Yield RNA transcript labeling kit (Enzo Diagnostics, 
NY, USA). Hybridization and scanning analysis were carried out as 
described before by Goda et al. (2002). 

Data analysis 

Image analysis and signal quantification were performed with 
QUANT array version 2.0 (GSI Lumonics, Oxnard, CA, USA). 
Background fluorescence was calculated on the basis of the fluor- 
escence signal of the negative control genes (the mouse nicotinic 
acetylcholine receptor epsilon-subunit (nAChRE) gene and the 
mouse glucocorticoid receptor homolog gene). A control template 
DNA fragment (TX803; Takara, Kyoto, Japan) was used as an 
external control to equalize hybridization signals generated from 
different samples. Gene-clustering analysis was performed with 
GENESPRING 5.0 (Silicon Genetics, San Carlos, CA, USA). 

A consensus motif search of promoter regions was performed at 
the Multiple Em for Motif Elicitation (http^/meme .sdsc.edu/me me/ 
website/meme.html). Analysis of consensus motifs was performed 
with DNASIS pro version 2.0 (Hitachi Software Engineering, 
Tokyo, Japan). Statistical significance was determined by x 2 test. 
The sources of DNA, sequences in promoter regions were used at 
the RIKEN Arabidopsis Genome Encyclopedia (http://rarge.gsc. 
riken.go.jp/). 



RNA gel blot analysis 

The RAFL cDNA probes whose sequences were determined by the 
Arabidopsis subspecies (SSP) sequencing consortium and RIKEN 
Genomic Sciences Center were amplified by polymerase chain 
reaction (PCR). RNA gel blot hybridization was performed as 
described before by Liu ef al. (1998). The RAFL cDNAs are available 
from RIKEN Bioresource Center. 



Acknowledgements 

We are grateful for the excellent technical support provided by 
Fumie Saito, Kyoko Yoshiwara, Chitose Kato, Ekuko Ohgawara, 
and Mie Yamamoto of Japan International Research Center for 
Agricultural Sciences (JIRCAS) Junko Ishida, Akiko Enju, Maiko 
Nakajima, Tetsuya Sakurai, Kei lida, and Masakazu Satou of 
Laboratory of Plant Molecular Biology, RIKEN and Kimihisa 
Ohshima of Saitama prefectural office. This work was supported 
in part by the Program for Promotion of Basic Research Activities 
for Innovative Biosciences (BRAIN). This work was also supported 



© Blackwell Publishing Ltd, The Plant Journal, (2004), 38, 982-993 



Cold-inducible downstream genes of DREB1A 993 



in part by a project grant from the Ministry of Agriculture, Forestry, 
and Fisheries, Japan. 

References 

Artschul, S.F., Gish, W., Miller, W., Myers, E.W. and Upman, DJ. 

(1990) Basic local alignment search tool, J. Mol. Biol. 215, 403- 
410. 

Bailey, T.L and Elkan, C. (1995) The value of prior knowledge in 
discovering motifs with MEME. Proc. Int. Conf. Intel!. Syst. Mol. 
Biol. 3, 21-29. . 

Baker, S.S., Wilhelm, K.S. and Thomashow, M.F. (1994) The 5'- 
region of Arabidopsis thaliana cor 15a has c/s-acting elements 
that confer cold-, drought- and ABA-regulated gene expression. 
Plant Mol. Biol. 24, 701-713. 

Bateman, A., Birney, E„ Durbin, R., Eddy, S.R., Finn. R.D. and 
Sonnhammer, E.L (1999) Pfam 3.1: 1313 multiple alignments 
and profile HMMs match the majority of proteins. Nucl. Acids 
Res. 27, 260-262. 

Bray, E.A. (1997) Plant responses to water deficit. Trends Plant Sci. 
2, 48-54. 

Chen, W., Provart, N.J., Glaze brook, J. et al. (2002) Expression 
profile matrix of Arabidopsis transcription factor genes suggests 
their putative functions in response to environmental stresses. 
Plant Cell, 14, 559-574. 

Fowler, S. and Thomashow, M.F. (2002) Arabidopsis transcrip- 
tome profiling indicates that multiple regulatory pathways are 
activated during cold acclimation in addition to the CBF cold 
response pathway. Plant Cell, 14, 1675-1690. 

Gilmour, SJ., Zarka, D.G., Stockinger, EJ., Salazar, M.P., 
Houghton, J.M. and Thomashow, M.F. (1998) Low temperature 
regulation of the Arabidopsis CBF family of AP2 transcriptional 
activators as an early step in cold-induced COR gene expression. 
Plant J. 16, 433^142. 

Gilmour, S.J., Sebort, A.M., Salazar, M.P., Everard, J.D. and Tho- 
mashow, M.F. (2000) Overexpression of the Arabidopsis CBF3 
transcriptional activator mimics multiple biochemical changes 
associated with cold acclimation. Plant Physiol. 124, 1854-1865. 

Goda, H., Shimada, Y., Asami, T., Fujioka, S. and Yoshida, S. (2002) 
Microarray analysis of brassinosteroid-regulated genes in 
Arabidopsis. Plant Physiol. 130, 1319-1334. 

Harmer, S.L, Hogenesch, J.B., Straume, M., Chang, H.S., Han, B., 
Zhu, T. f Wang, X., Kreps, J.A. and Kay, S.A. (2000) Orchestrated 
transcription of key pathways in Arabidopsis by the circa dian 
clock. Science, 290, 2110-2113. 

Ingram, J. and Bartels, D. (1996) The molecular basis of dehydra- 
tion tolerance in plants. Annu. Rev. Plant Physiol. Plant Mol. Biol. 
47,377-^03. 

Jaglo-Ottosen, K.R., Gilmour, S.J., Zarka, D.G., Schabenberger, O. 
and Thomashow, M.F. (1998) Arabidopsis CBF1 overexpression 
induces COR genes and enhances freezing tolerance. Science, 
280, 104-106. 

Jiang, C, lu, B. and Singh, J. (1996) Requirement of a CCGAC c/s- 
acting element for cold induction of the BN1 ?5gene from winter 
Brassica napus. Plant Mol. Biol. 30, 679-684. 

Kasuga, M., Liu, CL, Miura, S., Yamaguchi-Shinozaki, K. and 
Shinozaki, K. (1999) Improving plant drought, salt and freezing 
tolerance by gene transfer of a single stress-inducible transcrip- 
tion factor. Nat. Biotechnol. 17, 287-291. 

Liu, Q., Kasuga, M., Sakuma, Y., Abe, H., Miura, S., Yamaguchi- 
Shinozaki, K. and Shinozaki, K. (1998) Two transcription factors 



DREB1 and DREB2 with an EREBP/AP2 DNA binding domain 
separate two cellular signal transduction pathways in drought- 
and low-temperature-responsive gene expression respectively 
in Arabidopsis. Plant Cell, 10, 1391-1406. 

Logemann, E., Wu, S.C., Schroder, J., Schmelzer, E., Somssich, I.E. 
and Hahlbrock, K. (1995) Gene activation by UV light, fungal 
elicitor or fungal infection in Petroselinum crispum is correlated 
with repression of cell cycle-related genes. Plant J. 8, 865-876. 

Ohta, M., Matsui, K., Hiratsu, K., Shinshi, H. and Ohme-Takagi, M. 
(2001) Repression domains of class H ERF transcriptional repres- 
sors share an essential motif for active repression. Plant Cell, 13, 
1959-1968. 

Sakuma, Y., Liu, Q., Dubouzet, J.G., Abe, H., Shinozaki, K. and 
Yamaguchi-Shinozaki, K. (2002) DNA-binding specificity of the 
ERF/AP2 domain of Arabidopsis DREBs, transcription factors 
involved in dehydration- and cold-inducible gene expression. 
Biochem. Biophys. Res. Commun. 290, 998-1009. 

Schneider, T.D. and Stephens, R.M. (1990) Sequence logos: a 
new way to display consensus sequences. Nucl. Acids Res. 
18, 6097-6100. 

Seki, M., Narusaka, M., Abe, H., Kasuga, M., Yamaguchi-Shino- 
zaki, K., Carninici, P., Hayashizaki, Y. and Shinozaki, K. (2001) 
Monitoring the expression pattern of 1300 Arabidopsis genes 
under drought, and cold stress by using a full-length cDNA 
microarray. Plant Cell, 13, 61-72. 

Seki, M., Narusaka, M., Kamiya, A. etal. (2002) Functional annota- 
tion of a full-length Arabidopsis cDNA collection. Science, 296, 
141-145. 

Shinozaki, K. and Yamaguchi-Shinozaki, K. (1997) Gene expres- 
sion and signal transduction in water stress response. Plant 
Physiol. 115, 327-334. 

Shinozaki, K. and Yamaguchi-Shinozaki, K. (2000) Molecular 
responses to dehydration and low temperature: differences 
and cross-talk between two stress signaling pathways. Curr. 
Opin. Plant Biol. 3, 217-223. 

Shinwari, Z.K., Nakashima, K., Miura, S., Kasuga, M., Seki, M., 
Yamaguchi-Shinozaki, K. and Shinozaki, K. (1998) An Arabidop- 
sis gene family encoding DRE/CRT binding proteins involved 
in low-temperature-responsive gene expression. Biochem. 
Biophys. Res. Commun. 22, 161-170. 

Stockinger, EJ., Gilmour, SJ. and Thomashow, M.F. (1997) Ara- 
bidopsis thaliana CBF1 encodes an AP2 domain-containing 
transcriptional activator that binds to the C-repeat/DRE, a c/s- 
acting DNA regulatory element that stimulates transcription in 
response to low temperature and water deficit. Proc. Natl. Acad. 
Sci. USA, 94, 1035-1040. 

Tap, T„ Ohsumi, C, luchi, S., Seki, M., Kasuga, M., Kobayashi, M., 
Yamaguchi-Shinozaki, K. and Shinozaki, K. (2002) Important 
roles of drought- and cold-inducible genes for gaiactinot 
synthase in stress tolerance in Arabidopsis thaliana. Plant J. 
29, 417-426. 

Takatsuji, H. (1999) Zinc-finger proteins: the classical zinc finger 
emerges in contemporary plant science. Plant Mol. Biol. 39, 
1073-1078. 

Thomashow, M.F. (1999) Plant cold acclimation: freezing tolerance 
genes and regulatory mechanisms. Annu. Rev. Plant Physiol. 
Plant Mol. Biol. 50, 571-599. 

Yamaguchi-Shinozaki, K. and Shinozaki, K. (1994) A novel c/s- 
acting element in an Arabidopsis gene is involved in respon- 
siveness to drought, low-temperature or high-salt stress. Plant 
Cell, 6, 251-264. 



© Blackwell Publishing Ltd, The Plant Journal, (2004), 38, 982-993 



